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Abstract 

We present preliminary results of searches for charged B mesons decaying into the charmless three- 
body final states h^h^h^, where h = tt ov K, using 51.5 fb-i of data collected at the T{4S) 
resonance with the BABAR detector at the SLAC PEP-II asymmetric B Factory. No assumptions are 
made about intermediate resonances. We measure the branching fractions 13{B^ — > K^-k^tt^) = 
(59.2 ±4.7 ±4.9) x IQ-^ and B{B^ K^K^K^) = (34.7 ± 2.0 ± 1.8) x IQ-^, where the first error 
is statistical and the second error is systematic. In the same study, we do not observe significant 
signals for the final states B^ — > vr^vr^vr^ and B^ K^K^tt^, and therefore provide the 
confidence upper limits ^(5=^ 7r=^7r=F7r=^) < 15 x 10"^ and B{B^ K^K'^tt^) < 7 x lO"*'. 
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1 Introduction 



The study of charmless hadronic B decays is important to understand the phenomenon of CP 
violation in the Standard Model. There has been recent theoretical progress on using three-body 
decays to measure direct CP violation and to extract the Cabibbo-Kobayashi-Maskawa angle 7 Q. 
It is necessary to first observe these decays before such measurements can be made. We present 
updated preliminary results on the branching fractions of charged charmless three-body — > 
h^h^h^ decays, where h = it 01 K, with no assumptions about intermediate resonances and 
with open charm contributions subtracted. Charge conjugate initial and final states are assumed 
throughout this document, unless stated otherwise. 

2 The BaBAR detector and dataset 

The data used in this analysis were collected with the BABAR detector at the PEP-II asymmetric 
e~^e~ storage ring at SLAC. The data sample consists of 56.2 million BB pairs, corresponding to 
an integrated luminosity of 51.5 fb^^ collected at the T{AS) resonance (on-resonance) during the 
2000-2001 run. In addition, a total integrated luminosity of 6.4 fb~^ was taken at 40 MeV below 
the T(45) resonance (off-resonance), and was used to characterise the backgrounds from e~^e~ 
annihilation into light qq pairs. 

The BABAR detector is described in detail elsewhere ^J; the main parts relevant for the analysis 
of three charged particle final states are the tracking and particle identification sub-detectors. 

The 5-layer double-sided silicon vertex tracker (SVT) measures the impact parameters, angles, 
and transverse momenta of tracks down to 65 MeV/ c. Outside the SVT is a 40-layer drift chamber 
(DCH) that measures the transverse momenta of tracks from their curvature in the 1.5-T solenoidal 
magnetic field. The SVT and DCH also are used to determine the mean ionisation energy loss of 
tracks to help identify charged particles. The tracking system has a momentum resolution of 0.5% 
for a transverse momentum of 1.0 GeV/c. 

Surrounding the DCH is a detector of internally refiected Cherenkov radiation (DIRC), which 
provides charged hadron identification in the barrel region. Charged particles are identified by the 
Cherenkov angle Oc and the number of photons measured with the DIRC. The typical separation 
between pions and kaons varies from > 8 cr at 2.0 GeV/c to 2.5ct at 4.0 GeV/c, where a is the average 
resolution on 9c- The kaon selection efficiency is approximately 80%, which is the product of the 
particle identification algorithm efficiency with geometrical acceptance, for a pion mis-identification 
probability of 2%. 

The DIRC is surrounded by an electromagnetic calorimeter (EMC), made up of 6580 CsI(Tl) 
crystals, which is used to measure the energies and angular positions of photons and electrons with 
excellent resolution. The EMC is used to veto electrons in this analysis; the probability of mis- 
identifying electrons as pions is approximately 5%, while the probability of mis-identifying pions 
as electrons is below 0.3%. 

3 Analysis method 

The total branching fraction for each B^ h^h^h^ mode is measured over the whole Dalitz plot 
- all resonant and non-resonant contributions are included. A set of selection criteria is applied 
to reconstruct each mode separately. Each Dalitz plot is divided into many equal area cells to 
enable us to find the selection efficiency as a function of position in the Dalitz plot. We also 
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take into account continuum backgrounds and cross-feed between each signal mode from K and tt 
mis-identification. 

3.1 Candidate Selection 

We reconstruct B candidates from charged tracks, where each track must have at least 12 hits in the 
DCH, a maximum momentum of 10 GeV/c, a minimum transverse momentum of 100 MeV/c, and 
must originate from the beam-spot. We find three-charged track combinations to form the B can- 
didates, and we require that their energies and momenta satisfy kinematic constraints appropriate 
for B mesons. 

There are two variables we use for this, the first of which is the beam-energy substituted mass 
mES = \J {E"^ — p^)- The energy of the B candidate is defined as Ei, = {^s + Pq ■ pb)/Eo, where 
^/s and Eq are the total energies of the e~^e~ system in the centre-of-mass (CM) and laboratory 
frames, respectively, and po and pB are the momentum vectors in the laboratory frame of the e"^e~ 
system and the B candidate, respectively. The mss value should be close to the nominal B mass 
for signal events. 

The second variable we use is the energy difference between the reconstructed B candidate 
energy and the beam energy, AE = E'^ — \fs/2, where E*^ is the energy of the B candidate in the 
CM system. For this analysis, we assume the appropriate mass hypothesis for each charged track 
in a given decay mode under study in calculating /S.E. For signal events, AS should be centred at 
zero. The typical AE separation between modes that differ by substituting a kaon for a pion in 
the final state is 45 MeV. 

We use dE/dx information from the SVT and DCH, and the Cherenkov angle and number of 
photons measured by the DIRC for tracks with momenta above 700 McV/c, to identify charged 
pions and kaons. Kaons are selected with requirements made to the product of the likelihood 
ratios determined from these measurements. The likelihood ratio requirements are established by 
requiring the probability of mis- identifying pions as kaons be below 5%, up to a momentum of 
4.0 GeV/c. Pions are required to fail the kaon selection. We veto electron candidates by requiring 
that they fail a selection based on information from dE/dx, shower shapes in the EMC and the 
ratio of the shower energy and track momentum. 

Since we are only interested in charmless decays, we need to veto candidates that contain charm 
mesons. We remove B candidates when the invariant mass of the combination of any two of its 
daughter tracks (of opposite charge) is within 3 cr of the mass of , J /ijj or ijj{2S) mesons. Here, a 
is 10.0 MeV/c^ for and 15.0 MeV/c'^ for J/t/j and ip{2S). All possible kaon and pion combinations 
are tested for the veto, while only the K'^K~ and tt+tt" hypotheses are tested for the J/V' and 
ip{2S) vetoes, since the background from these decays is from leptonic decays, in which the leptons 
have been mis-identified as pions or kaons. The electron veto helps to reduce the combinatorial 
background from J/^p and ip{2S) decays that would otherwise pass the 3 a invariant mass veto. 

3.2 Background Suppression and Chciracterisation 

In addition to these candidate selection requirements, we need to suppress backgrounds from light 
quark and charm continuum production. We reduce these by imposing requirements on two topo- 
logical event shape variables computed in the T(4S') rest frame. 

The first event shape variable is the cosine of the angle Or between the thrust axis of the selected 
B candidate and the thrust axis of the rest of the event, i.e. all charged tracks and neutral particles 
not in the B meson candidate. For continuum backgrounds, the directions of the two axes tend 
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to be aligned because the daughters of the reconstructed candidate generally lie along the dijet 
axis of such events. Therefore, the distribution of |cos0t| is strongly peaked towards unity. For B 
events, the distribution of |cos0t| is isotropic because the decay products from the two B mesons 
are essentially independent of each other. The difference in the |cos0t| dependence allows us to 
discriminate between signal B decays and continuum background. 

The second event shape variable is a Fisher discriminant |^ JF, which is formed from the linear 
combination 

9 

J=' = ^aiXi (1) 

i=l 

of the input variables Xi. The available variables are the summed scalar momenta of all charged 
and neutral particles from the rest of the event within nine nested cones coaxial with the thrust 
axis of the B candidate. The coefficients Oi are chosen to maximise the separation between signal 
and background events. They are calculated for each signal mode separately using Monte Carlo 
signal and light quark continuum events. 

Figure [l| shows the Fisher distributions for T(4S') events, from the control sample B~ — > 
D^TT^jD^ K~TT~^ in Monte Carlo simulation and on-resonance data, and for background, from 
off-resonance data and light-quark continuum Monte Carlo events. 

The selection criteria for the event shape variables, shown in Table |l], is optimised separately 
for each signal mode to achieve maximum sensitivity for the branching fraction. 

Table 1: Selection requirements on the event shape variables for each signal mode. 



Signal Mode 


COS^T 


T 




< 0.575 


< -0.11 




< 0.700 


< -0.03 




< 0.725 


< 0.10 




< 0.875 


< 0.30 



Despite using the powerful event shape variables mentioned above, there are still significant 
backgrounds that must be explicitly subtracted to extract a signal. The residual background level is 
estimated from the observed number of events in a sideband region, located near to the signal region 
in the rriEs- AE plane, and then extrapolating into the signal region by using a multiplicative factor, 
R. We define R to be the ratio of the number of background candidates in the signal region to the 
number in the sideband region. In order to determine R, the shape of the background distribution 
as a function of tties is parameterised according to the ARGUS function Q , and is measured using 
the upper sideband in the AE variable in on-resonance data (0.1 < \AE\ < 0.25 GeV). A quadratic 
function is used to parameterise the background distribution as a function of AE. The ratio of the 
areas under the shape function in AE and rriES in the signal and sideband regions is equal to R. 
The uncertainty of the value of R is dominated by the uncertainty of the shape parameter for the 
ARGUS function. Off-resonance data give a consistent value of R. 

3.3 Branching Fraction Calculation 

As mentioned previously, the branching fractions for each signal mode are measured over the whole 
Dalitz plot, and each Dalitz plot is divided up into many cells so that the bin- by-bin variation of 
the selection efficiency can be found for each plot. 
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Figure 1: Normalised distributions of the Fisher discriminant for B~ D^'ir~ , K~'ir~^ Monte 

Carlo events (solid histogram), on-resonance D^it~ data (solid points), light quark continuum 
Monte Carlo events (dotted histogram) and off-resonance data (dotted points). 

Taking ej to be the Monte Carlo efficiency of reconstructing the signal in the i^^ bin in the 
Dalitz plot, the branching fraction for the signal mode is given by: 

where is the total number of BB pairs, R is the background extrapolation factor into the 

signal region, and Nu and N2i are the number of events observed in the signal and grand sideband 
(GSB) regions, respectively, for the Dalitz plot bin. No significant differences were found for 
the value of R in different regions of the Dalitz plot, so an average value is used for all bins. Si is 
the number of background subtracted signal events for the i Dalitz plot bin. 

The signal region is defined to be |mES - msl < 8.0 MeV/c^ and - (A^;) | < 60.0 MeV, 
where (Ai?) is the mean value of AE for on-resonance data for the calibration sample B~ — > 
D^-K~ — > K~'K~^, and uib is the nominal mass of the charged B meson The GSB region is 
defined to be 5.21 < mss < 5.25 GeV/c^ and \ AE - (AE) \ < 100.0 MeV. 

The probability of a kaon being mis-identified as a pion is 20%, which includes the efficiency 
of the particle identification algorithm and the geometrical acceptance. This means that there is 
significant cross-feed into the signal region from the decay mode that has one more kaon, which must 
be subtracted for each bin, i. To see how this is done, consider the branching fraction calculation 
for B^ K^K^K^. This channel has negligible cross-feed from the other modes, so we have 

B{B^ ^ K^K^K^) = ^ E - = f^' (3) 

BE i BB 
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where N^k is the total number of K^K^K^ events in the data. The channel K^K^tt^, on 

the other hand, has significant cross-feed from K^K^K^, and this must be subtracted from the 
signal. The branching fraction for this mode is given by 

B[B^ - K^K^^^) = ^ E ~ ^'""''^ = (4) 

^BB i ^BB 

where Nkktt is the total number of — > K^K^tt^ events in the data, and Si and refer 
to those quantities for B^ K^K^tt^. Here, e'l is the probability for reconstructing B^ 

events using the selection criteria for B^ K^K^n^. It is determined from Monte 
Carlo simulation by generating events uniformly in phase space and determining the cross-feed 
selection efficiency in each Dalitz plot bin. 

The branching fraction for B^ K^ir^n^ is given by 

B{B^ ^ K^n^TT^) = ^ E ~ ^^^-''^ - "^-^ = (5) 

^BB i ^BB 

where Nk-k-k is the total number of B^ — > K^tt^tt^ events in the data, and Si and refer 
to the number of background subtracted events in the signal region and selection efficiency for 
B^ — > K^TT^TT^, respectively. Here, e-' is the probability for reconstructing B^ — > K^K^n^ 
events using the selection criteria for B^ — > K^tt^tt^. This channel has some contamination 
from candidates falling outside the 3 a invariant mass veto, which must be subtracted. This is 
represented by the term n^ij, which is the number of events that is expected to populate the 

bin in the Dalitz plot. This is estimated by finding the probability of reconstructing Monte 
Carlo events, for each bin i, using the selection criteria for B^ —>■ K^tt^tt^, and multiplying this 
by the measured branching fraction for the mode and the total number of BB pairs, N^^. 
The total expected number of B~ — > D^tt~,D^ K~Tr~^ events that must be subtracted across 
the full Dalitz plot for this channel is 47 it 8. The values of noi are non-zero only for bins close to 
the resonance bands. 

Finally, the branching fraction for B^ tt^vt^tt^ is given by 

^BBi ^BB ' 

where N^.,^ is the total number of B^ — > Tr^vr^vr^ events in the data, e" is the probability for 
reconstructing B^ — > K^tt^tt^ events using the selection criteria for B^ — > tt^tt^tt^, and Si and 
ei refer to the number of background subtracted events in the signal region and selection efficiency 
for B^ — > vr^vr^vr^, respectively. Again, there are some B- D^7r-,D^ K-7r+ events (23 ± 5) 
that pass the selection criteria for this channel, because the kaon from the decay is mis-identified 
as a pion, and where the invariant mass of the meson lies outside the 3 a invariant mass window. 
This background is subtracted from the signal. 

In addition to the cross-feed where only one of the kaon tracks is mis-identified as a pion, there 
can also be cross-feed where either two kaons are mis-identified as pions (probability of 4%) or when 
one of the pions is mis-identified as a kaon (probability of 2%). These are smaller, second-order 
cross-feed effects, and so it is adequate to subtract the average number of events over the whole 
Dalitz plot. If Ux is the average number of second-order cross-feed events that has to be subtracted 
(i.e. the number of events reconstructed divided by the appropriate cross- feed efficiency), then we 
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finally have 

where Nx is the total number of events from the channel that contributes to most of the (first-order) 
cross-feed, e.g. A^^^ = N^^^^ for — > TT^ir^ir^. 



Table 2: Efficiencies and cross-contamination probabilities between the signal modes derived from 
Monte Carlo samples. For example, the probability that an event K^tt^tt^ will be reconstructed 
as 7r±7r=F7r=^ is (1.7 ± 0.1) x 10"^. 



Selection Criteria 
Hypothesis 


Input Decay Mode 












(15.3 ±0.2) X 10^^ 
(3.6 ±0.4) X 10-3 
(0.0 ±0.2) X 10-3 
(0.0 ±0.2) X 10-3 


(1.7±0.1) X 10-^ 
(15.1 ±0.2) X 10-2 
(2.9 ±0.4) X 10-3 
(0.0 ±0.2) X 10-3 


(1.4 ±0.9) X lO"'' 
(3.2 ±0.2) X 10-2 
(17.7 ±0.3) X 10-2 
(1.7 ±0.2) X 10-3 


(1.1 ±3.2) X 10-^ 
(4.0 ± 1.7) X 10-"* 
(5.5 ±0.2) X 10-2 
(21.6 ±0.3) X 10-2 



The Dalitz plot for each signal mode is divided into cells with equal area 1.0 GeV^, and large 
samples of Monte Carlo signal events are used to obtain the signal and cross-feed selection efficiencies 
across each Dalitz plot. Table |^ shows the signal and cross-feed selection efficiencies for the modes, 
averaged over the Dalitz plots. 



4 Experimental Uncertainties 

There are several sources of uncertainty for the branching fraction measurements, which come from 
the various terms in Eq. |^. The statistical uncertainties come from the number of events observed 
in the signal and GSB regions, A'^i and N2- The factor R, found independently for each mode, 
has a systematic uncertainty arising from the uncertainty in the fitted ARGUS shape parameter, 
^. The main sources of S-related background are decays (the n^j term) and cross- feed from 
the other signal modes (the Nx and Ux terms), owing to kaon and pion mis-identification. We deal 
with these by explicitly subtracting them from the signal. The uncertainty in the number of 
events that have to be subtracted comes from the uncertainty on the published measured branching 
fractions the number of BB events (mentioned below) and the selection efficiency. 

Since there are a lot of terms used to calculate the branching fraction, it is worthwhile to go 
through what uncertainties they contribute to the end result. If we let Xi represent the term within 
parenthesis in Eq. ^, then the fractional uncertainty on the branching fraction is 



(8) 



where Se is the fractional systematic error for the efficiency that comes from differences between 
Monte Carlo simulation and on-resonance data, shown in Table |3|. Going through the terms for X^, 
we have 
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+ 2^|^Aei -2 i J +(Ana;)^ (9) 

The first two terms on the right hand side (R.H.S.) are the statistical uncertainties on the number 
of events in the signal and GSB regions, while the third term represents the systematic variation 
for the background extrapolation factor, R. The uncertainties from the cross-feed subtraction 
are represented by the next three terms, while the penultimate term is that for the bin-by-bin 
uncertainty for the efficiency. The last term is the uncertainty for the number of second-order 
cross-feed events. The various sources of error mentioned above will be shown in detail for each 
branching fraction result. 

The uncertainties on the signal efficiencies and cross-feed probabilities are the combination of 
statistical errors on the number of events selected in the Monte Carlo samples relative to the total 
number generated, as well as systematic uncertainties arising from the difference between Monte 
Carlo simulation and on-resonance data. 

The average fractional Monte Carlo statistical uncertainties of the signal efficiencies per Dalitz 
plot bin (Aej/eO are 7.0% for 7r±7r^7r±, 9.1% for K^vt^tt^, 9.5% for K^K^t^^ and 7.4% for 

The Monte Carlo simulation is subject to systematic uncertainties from tracking and particle 
identification efficiencies. Residual differences in the tracking selection efficiencies between on- 
resonance data and Monte Carlo simulation contributes a fractional uncertainty of 0.8% on the 
efficiency per track. This uncertainty is added coherently for all three tracks used to reconstruct 
each B meson. Both the electron veto and kaon selections have fractional systematic uncertainties 
of 1.0%. These uncertainties are added coherently. 

Possible differences between the behaviour of Monte Carlo simulation and on-resonance data 
are also examined for the Fisher distributions used to discriminate signal B decays from light quark 
continuum events. The control samples B^ D^h^,D^ h~h'^, where /i = vr and/or K, which 
have similar kinematics to the signal modes, are used to compare the signal Fisher distributions 
between on-resonance and Monte Carlo data, using the Fisher coefficients derived from Monte Carlo 
signal and qq samples. The choice for h is made such that the final state of the control sample 
decay has the same number of kaons and pions as those for the signal mode. 

The Fisher distribution for off-resonance data agrees with that for light quark continuum Monte 
Carlo events, which can be seen in Fig. ^ There are very slight differences between the Fisher 
distributions for Monte Carlo signal events and on-resonance data. To quantify this difference, 
the Fisher distributions are fitted to double Gaussian functions. The differences in the mean and 
width values between on-resonance and Monte Carlo D^h events are used to shift and scale the 
Fisher distributions for the signal Monte Carlo modes. The change in the selection efficiency gives 
an estimate of the correction factor necessary for the requirement on the Fisher variable used in 
the selection for each mode, which is found to be very small (approximately 1%). The systematic 
uncertainty on this correction is found by varying the parameters of the Fisher distribution for the 
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signal Mode Carlo modes by the (scaled) uncertainties found for the fitted parameters for the 
control sample. 

The resolutions for AE and mss in data differ by a negligible amount from Monte Carlo 
predictions. The main source of uncertainty arises from a +7 MeV shift in the mean value of AE 
observed in the control sample B~ — > D^7r~,D^ K~Tr~^ in on-resonance data, which we correct 
for in the Monte Carlo. This contributes a fractional systematic uncertainty of 1%. 

Table ^ gives a summary of the systematic uncertainties to the efficiency for each mode (not 
including the Dalitz plot variation). The fractional uncertainties for the Dalitz plot variation for 
the cross-feed probabilities {Ae'- /e") are approximately 30%. 

Table 3: Fractional systematic uncertainties for the Monte Carlo efficiency for each signal mode. 
The uncertainties are added in quadrature in the total. 



Source of Uncertainty 


Fractional error on efficiency (%) 










Tracking 


2.4 


2.4 


2.4 


2.4 


Fisher Discriminant 


2.1 


0.7 


0.5 


1.7 


Particle Identification 


6.0 


5.0 


4.0 


3.0 


AE and tties 


1.0 


1.0 


1.0 


1.0 


Total 


6.8 


5.6 


4.7 


4.2 



Finally, there is a systematic uncertainty on the overall normalisation, N^^, which is obtained 
from a dedicated study to find the number of B mesons produced in the data sample. This is found 
to have a systematic uncertainty of 1.5%. 

5 Physics results 

Our preliminary measurements of the branching fractions for the signal modes are summarised in 
Table ^. The top few rows of this table show the total number of events observed in the signal and 
GSB regions, as well as the average signal reconstruction efficiencies for each mode and the values 
of the background extrapolation factor R. 

The row labelled 1) shows the sum over Dalitz plot bins of the number of events observed in 
the signal region divided by the signal efficiency. The error on these quantities is the first error 
shown in Eq. and only includes the uncertainty in the number of signal events, Nu. 

The next row, labelled 2), shows the sum over Dalitz plot bins of the expected number of 
background events divided by the signal efficiency. The errors shown for these values are the 
second and third terms on the R.H.S. of Eq. |9|, respectively. They correspond to the statistical 
uncertainty in N2i, and the systematic error for R, which is dominant. 

Row 3) shows the expected number of cross- feed events (from K/ir mis- identification) . The 
errors on these quantities represent the fourth and fifth terms on the R.H.S. of Eq. ^ only. Note 
that the K^K^K^ mode does not have a cross-feed term, since this is negligible. 

The expected number of events passing the selection criteria for B^ —>■ n^jr^Tr^ and B^ 
K^TT^TT^ is shown in row 4), where the error for each value is the sixth term on the R.H.S. of 
Eq. |. 

The second-order cross- feed terms, shown in row 5). There are only entries for B — s- 

K^tt'^tt^, from B^ K^K^K^ events, and for B^ K^K^tt^, from B^ K^tt'^tt^ cross- 
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feed. The errors for these values are dominated by the uncertainties in the second-order cross- feed 
probabihties. Note that the Ux term for — > K^tt^-k^ is negative, which compensates for the 
extra cross-feed background of — > K^K^K^ events that is mis-identified as B^ K^K^ir^, 
and then in turn passes the selection criteria for B^ — > K^tt^tt^. 

The various contributions to the signal and background terms are shown in row 6). The first 
uncertainties are the combination of the statistical errors for the number of events in the signal and 
GSB regions - they correspond to the sum in quadrature of the error for row 1), and the first error 
in row 2). The second error for the entries in row 6) corresponds to the quadrature sum of all the 
other systematic errors from rows 2) to 5). The third error for row 6) is that from the penultimate 
term on the R.H.S. of Eq. ^, i.e. the uncertainty for the selection efficiency for each of the Dalitz 
plot bins separately (not the average). The last error for row 6) is just the fractional systematic 
uncertainty for the efficiency correction factors, given in Table ^ 

The last row in Table ^ shows the branching fraction results, where the first uncertainties are the 
statistical errors on the number of signal and background events, while the second uncertainties 
are the sum in quadrature of all systematic errors. The dominant systematic uncertainty for 
B^ K^K^K^ is the systematic correction factor for the efficiency, while for B^ — > vr^vr^vr^ 
and B^ K^K^tt^ , the background extrapolation factor R dominates. For B^ — > K^tt^tt^, 
both uncertainties contribute equally to the systematic error. 

As a consistency check, the branching fraction for the control sample B~ D^tt^ , K~tt~^ 
is measured to be (180 it 4 it 11) x 10^^, which agrees with the previously measured value of 
(203 ±20) X 10-6 @. 

Figures ^ to ^ show the AE and ttt-es distributions for the signal region for each of the modes. 
Each plot shows the expected levels continuum and BB background (solid and dashed lines, re- 
spectively). 

Figures ^ to ^ show the unbinned Dalitz plots for the signal modes in the GSB and signal 
regions, where no efficiency corrections have been applied. Only the upper half of the symmetrical 
Dalitz plot is shown for the B^ vr^vr^vr^ and B^ — > K^K^K^ channels, where the x and y 
axes show the minimum and maximum values of the Dalitz plot variables, respectively. 

There are clear signals for the modes B^ K^ir^n^ and B^ — > K^K^K^. No signal is 
observed for B^ — > K^K^tt^ , and the result for B^ — > vr^vr^vr^ is interpreted as an upper limit 
on the branching fraction, although there is a positive excess of signal events with 2.2 a significance. 
Since there are a large number of events in the selected samples, we can assume that the number 
of signal and background events observed in the signal region are Gaussian distributed. The 90% 
C.L. upper limits are computed using the standard prescription for a one-sided confidence interval 
from a Gaussian distributed measurement, i.e., 

BuL = l3 + 1.28A^, (10) 

where B is the estimated branching ratio and AB is its standard deviation. Here, however, we take 
AB to be the total error (quadratic sum of statistical and systematic uncertainties). 
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Table 4: Branching fraction results for on-resonance data. The uncertainties for each term are explained in the text. 



Signal Mode 










No. of events in signal region, J^i 


951 


1269 


573 


603 


No. of events in GSB, J^i ^21 


5470 


4652 


3239 


1100 


Average signal efficiency (%) 


15.3 ± 1.1 


15.4 ±0.9 


18.3 ±0.9 


22.5 ± 1.0 


Background factor, R 


0.145 ±0.006 


0.153 ±0.006 


0.150 ± 0.006 


0.159 ±0.010 


1) EiA^uAi 


5839 ± 212 


8055 ± 255 


3414 ± 156 


2734 ± 111 


2) EiRN2i/ei 


4812 d= 73 ± 193 


4434 ± 73 ± 171 


2802 ± 54 ± 111 


780 ± 23 ± 47 


3) EiN.e'^/ei 


391 ± 8 ± 2 


14 ± 1 ± 1 


435 ± 5 ± 7 




4) No. of D*^ events, Yl,i''^Dil (■% 


157 ±27 


401 ± 50 






5) 2^^-order cross-feed, Ux 




-122 ± 54 


57 ±11 




^ {Nu-RN2i-N^e'^-nDi) 
6) l^i ^. -rix 


478 ± 224 ± 195 
±34 ± 33 


3328 ± 266 ± 186 
±56 ± 186 


120 ± 166 ± 112 
±22 ±6 


1954 ± 114 ± 47 
±13 ± 82 


Branching Fraction (xlO~^) 
Statistical Significance (cr) 
90% Upper Limit (xlO"^) 


8.5 ±4.0 ±3.6 
2.2 
< 15 


59.2 ±4.7 ±4.9 
> 6 


2.1 ±2.9 ±2.0 
0.9 
< 7 


34.7 ± 2.0 ± 1.8 
> 6 




Figure 2: On-resonance signal region AE (a) and tit-es (b) distributions for B — > vr tt^tt . 
The solid lines show the expected level of continuum background, using appropriately normalised 
background shapes from the sideband regions in on-resonance data. The dotted lines show the 
expected level of BB background, which is obtained from the sum of Gaussian distributions from 
Monte Carlo estimated cross-feed and D^ir events, each normalised to the number of events observed 
in on-resonance data that passed the selection criteria for B^ vr^vr^Tr^. 




Figure 3: On-resonance signal region AE (a) and ttt-es (b) distributions for B^ K^tt^tt^. 
The solid lines show the expected level of continuum background, using appropriately normalised 
background shapes from the sideband regions in on-resonance data. The dotted lines show the 
expected level of BB background, which is obtained from the sum of Gaussian distributions from 
Monte Carlo estimated cross-feed and D^ir events, each normalised to the number of events observed 
in on-resonance data that passed the selection criteria for B^ K^tt^tt^. 
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Figure 4: On-resonance signal region AE (a) and mss (b) distributions for K^K^tt^. 
The solid lines show the expected level of continuum background, using appropriately normalised 
background shapes from the sideband regions in on-resonance data. The dotted lines show the 
expected level of BB background, which is obtained from the sum of Gaussian distributions from 
Monte Carlo estimated cross-feed events, each normalised to the number of events observed in 
on-resonance data that passed the selection criteria for B^ —>■ K^K^tt^. 
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Figure 5: On-resonance signal region AE (a) and ttt-es (b) distributions for B^ K^K^K^. 
The solid lines show the expected level of continuum background, using appropriately normalised 
background shapes from the sideband regions in on-resonance data. 
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Figure 6: Unbinned Dalitz plots for on-resonance data for vr^vr^vr^ for GSB region (a) and 

signal region (b). No efficiency corrections have been applied, and the open charm contributions 
are included in the plots. 
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Figure 7: Unbinned Dalitz plots for on-resonance data for K^tt^tt^ for GSB region (a) and 

signal region (b). No efficiency corrections have been applied, and the open charm contributions 
are included in the plots. 
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Figure 8: Unbinned Dalitz plots for on-resonance data for K^K^tt^ for GSB region (a) and 

signal region (b). No efficiency corrections have been applied, and the open charm contributions 
are included in the plots. 



6 Summary 

We have obtained preliminary branching fractions for — > K^tt^tt^ and B^ K^K^K^ over 
the whole Dalitz plot, and have determined conservative 90% upper limits for B^ vr^Tr^vr^ and 
B^ K^K^TT^. The results are summarised in Table |5|, where the results from BELLE are 
also included for comparison. 



Table 5: Branching fraction results from BABAR and BELLE. 



Decay mode 


BABAR 


BELLE 




< 15 X 10"*^ 






(59.2 ±4.7 ±4.9) x 10"^ 


(58.5 ±7.1 ±8.8) X 10-'^ 




< 7 X 10^« 


< 21 X 10^^ 




(34.7 ±2.0 ±1.8) X lO"*^ 


(37.0 ±3.9 ±4.4) X lO"*^ 
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Figure 9: Unbinned Dalitz plots for on-resonance data for K^K^K^ for GSB region (a) and 

signal region (b). No efficiency corrections have been applied, and the open charm contributions 
are included in the plots. 
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